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UDC  629.12.053:681.2.082.16 
Gyroscopic  Control  of  Ship  Heading 

907F0172A  Leningrad  SLTOSTROYENIYE  in  Russian  No  9,  Sep  89  pp  25->26 
[Article  by  A.  Yu.  Panov  and  Ye.  I.  Shapkin] 

[Text]  The  traditional  devices  for  controlling  the  heading  of  a 
ship — rudders  and  rotary  nozzles — have  a  number  of  deficiencies,  the 
main  ones  of  which  are: 

poor  controllability  and  loss  of  it  at  low  speeds,  with  no  running 
and  in  reverse,  which  results  in  an  increase  of  the  mass  and  overall 
dimensions  of  the  rudders  and  results  in  installation  of  steering 
devices ; 

large  power  losses  of  the  ship  power  plant  when  ships  are  moving 
with  transposed  rudders  (rotary  nozzles)  due  to  the  occurrence  of  a 
longitudinal  component  of  the  hydrodynamic  force  that  slows  the  ship, 
and  up  to  40  percent  or  more  of  speed  is  lost  in  turning  [1,  2]; 

frequent  breaks  of  the  rudders  and  rotary  nozzles  when  moving 
through  clogged  and  shallow  channels  and  under  icing  conditions,  which 
results  in  large  losses  on  repair,  related  to  placing  the  ship  in 
drydock  or  a  caisson. 

The  indicated  deficiencies  are  corrected  by  using  gyroscopic  devices  for 
controlling  the  heading  of  the  ship.  It  is  obvious  from  a  patent  survey 
of  the  problem  that  studies  are  being  conducted  in  Great  Britain  [3], 
Japan  [4],  and  the  USSR  [5]  on  gyroscopic  ship  control  devices.  The 
physical  basis  of  the  method  of  controlling  the  heading  is  the 
capability  of  a  one-axis  gyroscope  of  generating  a  g^^roscopic  moment 
with  respect  to  the  hull  of  the  ship.  Let  a  one-axis  gyroscope  1 

(Figure  1)  rotate  about  the  principal  axis  at  angular  velocity  y,  the 
principal  axis  itself,  strengthened  in  a  gimbal  2,  rotates  at  angular 

— f 

velocity  u/\  .  According  to  the  law  of  mechanics  of  the  moment  of 
momentum,  we  have: 
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dK(y  -►  -► 
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(1) 


where  Ko  is  the  moment  of  momentum  of  a  gyroscope  with  respect  to  the 
center  0,  U  is  the  velocity  of  the  end  of  the  vector  of  the  moment  of 

-4  -i 

momentum,  mo(F®)  is  the  moment  of  external  forces  F®,  acting  on  the 
gyroscope  through  the  bearings  3. 


Figure  1.  Power  Diagram  of  Gyroscope 

-4  -4 

The  gyroscope  in  turn  acts  on  the  bearing  seats  3  with  forces  Fj  and  Fo , 

which  form  a  pair  with  moment  m^p .  This  moment  is  transmitted  to  the 
ship  hull,  vdiich  also  causes  it  to  turn  along  the  heading. 

Thus,  to  change  the  ship’s  heading,  the  principal  axis  of  the  gyroscope 
must  be  rotated  about  the  axis  parallel  to  the  main  plane  of  the  ship 
[5],  The  resulting  gyroscopic  moment  is  calculated  by  the  formula 
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(2) 


where  I  is  the  moment  of  inertia  of  the  gyroscope  with  respect  to  the 
principle  axis,  m  is  the  mass  of  the  gyroscope,  R  is  its  radius  and  y  is 
the  angle  between  the  principal  axis  of  the  gyroscope  and  the  base 
plane . 

An  experiment  was  conducted  on  the  model  of  a  ship  to  check  the 
theoretical  proposition.  The  specifications  of  the  model  of  a  ship  with 
gyroscope  are  as  follows;  length  of  the  hull  of  =  1.135  m,  width  of 

the  hull  of  B  =  0.196  m,  draft  of  T,  =  0.052  m,  mass  of  m  =  6.45  kg, 

mass  of  the  gyroscope  of  m^.  =  0.2  kg,  moment  of  momentum  of  Kc  =  0.392 

kgm2/s,  and  rotational  frequency  of  n^,  =  2,250  rpm.  The  dependence  of 

the  angular  rate  of  turn  iiiz  and  course  angle  of  the  model  on  the 
angular  rate  of  turn  (nidder  angle)  of  the  principal  axis  of  the 
gyroscope  lyi  (Figure  2)  was  found  as  a  result  of  the  tests. 


Figure  2.  Dependence  of  Angular  Velocity  of  Ship  u>z 
and  of  Course  Angle  f  on  Angular  Rate  of  Transfer  of 
Principle  Axis  of  GjToscope  ij\ 


KEY: 

1.  Radians  2.  Radians/s 

The  technical  capability  of  mounting  a  heading  control  gyroscope  on 
various  cleisses  of  ships  was  determined  from  the  condition  of  the 
equality  of  moments  that  turn  the  ship  on  the  heading,  vAiich  are  created 
by  the  control  devices  and  gyroscopes  on  ships,  resipectively .  The 
hydrodynamic  forces  and  moments  occurring  on  the  rudders  and  rotary 
nozzles  were  calculated  by  known  methods  [1,  6]. 
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The  gyroscope  is  selected  by  the  gyroscopic  moment  m  ,  which  guarantees 

Bp 

the  required  controllability  to  the  ship.  The  rotational  frequency  of 
the  gyroscope  and  its  radius  are  determined  by  the  stresses  in  the  rim 
of  the  flywheel 


C  s*  pv2^ 


(3) 


ubere  v  is  the  circumferential  velocity  on  the  rim  of  the  flywheel,  m/s, 
and  p  is  the  density  of  the  flywheel  material,  kg/m3. 

The  stresses  (r  are  limited  by  the  yield  point  ir^  of  the  flywheel 

material .  40KCrNiMo  alloy  is  used  as  the  material  [ 7 ] .  The  use  of 
ccmiposite  flytdieels  [8]  guarantees  a  reduction  of  the  overall  dimensions 
and  weight  characteristics  of  gjToscopes. 

The  power  of  the  electric  motor  that  rotates  the  gyroscope  is  determined 
from  calculation  of  the  frictional  losses  of  the  fljTdieel  against  air 
(ventilation  losses)  and  in  their  bearings.  The  ventilation  losses  (kW) 
for  open  flyv4ieels  are  determined  by  the  empirical  formula  [8] : 


A’b  =  (5,7£)  +  291)  D'/!3.  10-3, 


(4) 


where  D  is  the  diameter  of  the  flywheel,  m,  L  is  the  width  of  the 
flywheel,  m,  and  n  is  the  rotational  frequency,  rpm. 

These  losses  are  considerable  for  large-diameter  flyvdieels  (they  are 
comparable  to  the  power  of  the  steering  devices).  However,  there  are 
effective  methods  of  reducing  ventilation  losses:  installation  of  the 
flyidieel  in  a  closed  housing,  creation  of  a  vacuum  in  it,  or  filling  it 
with  light  gases  [9]. 

The  gyroscope  is  controlled  by  an  electromechanical  drive,  consisting  of 
an  electric  motor  and  reduction  gear.  Controlling  the  electric  motor, 
the  navigator  thus  controls  the  ship.  The  power  of  the  electric  motor 
of  the  drive  (kW)  is  determined  by  the  formula 


A^„  =  Af-fir;, 


(5) 


where  7  =  u>i  is  the  angular  rotational  velocity  of  the  principal  axis  of 
the  gyroscope,  rad/s,  ij  is  the  efficiency  (kN)  of  the  reduction  gear,  M 
is  the  moment  on  the  rotational  axis  of  the  gyroscope  gimbal ,  determined 
by  the  approximate  formula: 
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Design  Parameters  of  Gyroscopes  for  Different  Classes  of  Ships 
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KEY: 
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M  =  cos  7.  ( 0 ) 


vrfiere  Iz  is  the  mcanent  of  inertia  of  the  g>Toscope  with  respect  to  the 

principal  axis,  f  is  the  angular  rate  of  turn  of  the  ship,  and  0  is  the 
course  angle. 

The  precise  formula  for  determination  of  moment  M  is  presented  in  [10]. 
Calculations  by  formulas  (4)-(6)  for  an  open  fljTdieel  D  =  0.75  m  in 
dimeter  and  L  =  0.46  m  wide  at  rotational  frequency  of  n  =  4,500  rpm, 
v^ich  determine  the  controllability  of  a  motor  ship  18  m  long,  indicate 
values  of  N-  a  5  kW  and  N„  2:  2  kW. 

For  ships  operating  in  heavy  seas,  two-flyv4ieel  gyroscopes  must  be  used 
(the  calculated  mass  of  the  gyroscope  is  redistributed  between  the  two 
flyvdieels ) .  The  vectors  of  the  moments  of  momentum  of  the  gyroscopes  in 
the  initial  ^sition  are  established  in  opposite  directions,  while 
linear  velocities  of  the  same  direction  are  assigned  to  their  ends. 

Thus,  the  harmful  lateral  gyroscope  moments  that  occur  upon  rolling  of 
the  ship  i^d  that  have  a  negative  influence  on  its  stability  are 
compensated . 

The  results  of  calculating  the  structural  parameters  of  control 
gyroscopes  are  presented  in  the  table.  The  range  of  application  of 
control  gyroscopes  is  obvious  from  the  calculations:  this  is  a  ship  of 


Class  of  ship 
Displacement ,  t 
Speed,  km/hr 
Calculated  length,  m 
Draft,  m 

Mcanent  on  ship  hull,  created  by  trauiitional  control  devices, 
kN*m 

Parameters  of  gyroscopes  that  create  same  moment  on  ship 
hull 

Rotational  frequency  of  gyroscope,  rad/s 
Radius  of  gyroscope,  m 
Mass  of  gyroscope,  t 

Mass  of  gyroscope  in  percent  of  displacement 
Dry-cargo  motor  ship  of  "Yakutsk"  class 
Dry-cargo  motor  ship  of  "Okskiy-50"  class 
Dry-cjargo  motor  ship  of  "Oka-l"  class 
Tanker  of  "Samotlor"  class 

Self-propelled  dirt-hauling  barge  with  capacity  of  300  m^ 
Pusher  tug  of  OT-2101  class 
Tug  of  "Ulan-Ude"  class 
Harbor  tug  of  BT-861  class 
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displacement  up  to  2,000  t,  operating  on  rivers  with  complex  channel. 

The  use  of  these  devices  on  ships  of  the  port  and  commercial  fleet  is 
especially  advantageous.  The  masses  of  gyroscopes  for  ships  of  various 
classes  comprise  an  average  of  0.3-0. 7  percent  of  the  mass  of  the  loaded 
ship,  vrfiereas  this  indicator  reaches  0. 7-1.0  percent  for  rudders.  It  is 
expected  that  the  mass  of  control  gyroscopes  with  regard  to  the  mass  of 
the  gimbal,  substructures  and  electric  motors  will  be  comparable  to  the 
mass  of  existing  ship  control  devices. 

In  conclusion,  we  note  the  main  advantages  of  ship  heading  control 
gyroscopes : 

they  guarantee  the  required  controllability  at  low  speeds,  in 
reverse  and  in  the  absence  of  running  speed,  where  other  control  devices 
are  hardly  effective; 

the  absence  of  interaction  with  the  water  guarantees  high 
efficiency  of  these  devices,  conserv^ation  of  fuel,  and  an  increase  of 
running  speed; 

location  inside  the  ship  hull  guarantees  higher  reliability  and 
rei^airabili  ty ; 

ships  equipped  with  gyroscopic  control  are  stable  in  heading,  since 
the  gyroscope  spontaneously  processes  and  generates  a  moment  that  holds 
the  ship  on  course  when  the  ship  is  affected  by  perturbing  forces. 
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UDC  5S1 .7.082.5 

Measurement  of  the  Amplitudes  and  Phases  of  Surface  Acoustic  Waves  Using 
an  Automated  Laser  Installation 

907F0217  Ashkhabad  IZVF.STTYA  AKADF.MTI  NAUK  TURKMENSKOY  SSR:  SF.RTYA 
FIZIKO-TFKHNICHESKIKR,  KHIMTCHESfCIKH  I  GEOI.OGfCHESKIKH  NAUK  in  Russian 
Mo  6,  Nov-Dee  89  pp  79-81 

[.Article  by  S.  S.  Sukhanov,  G.  Garyagdyyev,  K.  I.  Britsyn,  V.  G.  Grits,  K. 
Nurmukhamedov,  V.  V.  Semenov,  C.  Ch.  Sushinski] 

[Text] 

Optical  sensing  methods  in  current  u.se  make  it  possible  to  measure  the  amplitude 
and  phase  of  surface  acoustic  waves  (SAWs)  in  various  acoustico-electronic  devices 
without  disturbing  the  medium  under  study.  The  devices  studied  include  delay 
lines,  band  and  matched  filters,  correlators,  resonators,  etc.  One  needs  to  know  the 
parameters  of  SAWs  to  monitor  physical  processes  occurring  in  the  material  of 
these  products  and  to  improve  the  characteristics  of  the  products  and  instruments 
which  are  being  developed. 

Known  installations  for  measuring  the  parameters  of  SAWs  are  bulky  and  difficult 
to  tune.  They  are  not  adapted  to  work  in  a  wide  band  of  frequencies  or  to  monitor 
the  parameters  of  products  with  S.AWs  during  their  production  [1]. 

The  automated  laser  installation  described  below  measures  the  amplitude  and  phase 
profiles  of  SAWs.  It  is  to  a  great  extent  devoid  of  the  drawbacks  of  known 
installations  of  this  type.  Measures  taken  in  the  circuitry  and  construction  of  this 
device  have  made  it  possible  to  obtain  a  high  sensitivity  to  small  SAW  amplitudes, 
measurement  stability,  simplicity  of  u.se,  and  the  capacity  to  work  in  a  wide  range  of 
SAW  frequencies. 

The  optical  schematic  of  the  installation  is  shown  in  Eigure  1.  The  optical  quantum 
generator  1,  type  LG-38  through  a  prism  with  two  reflecting  surfaces  2  is  equipped 
with  a  heterodyne  channel  of  the  installation  with  monochromatic  radiation  (\o  = 
0.63  /tm);  in  acoustico-optical  modulator  3  (a  Dragg  cell)  there  is  a  spatial  and 
frequency  separation  of  the  laser  radiation  by  orders  of  Bragg  diffraction.  Two 
orders  of  diffraction  are  used,  the  zero  and  the  first.  The  forward  ray  of  the  zero 
order  is  focused  by  micro-objective  5  onto  the  surface  of  the  product  under  study 
6  with  SAWs.  The  reflected  ray,  .which  carries  information  on  the  amplitude  and 
phase  of  the  SAWs,  returns  to  the  acoustico-optical  modulator,  and  diffracting  in  it, 
falls  on  the  receiving  window  of  the  photodetector  9. 

The  ray  of  the  quantum  generator  of  the  order  of  diffraction  is  focused  by  the 
micro-objective  8  onto  the  surface  of  the  mirror  7  and,  returning  through  the 


9 


Figure 

sound. 


The  optical  schematic  of  the  installation,  a.  photoreceiver;  b.  node;  c. 


acoustico-optical  modulator,  also  enters  the  receiving  window  of  the  photodetector, 
forming  the  heterodyne  channel  of  the  installation. 

Tn  the  photodetector,  which  is  operating  in  nonlinear  mode,  there  is  a  transformation 
of  optical  signals  into  radio  frequency  signals: 

2/b.  2/s  -  /c.  2/b  +  /^,  (1) 


CA /c  ’S  the  frequency  of  the  acoustic  excitation  of  the  Bragg  cell;  is  the 
SAW  frequency.  These  signals  carry  information  on  the  amplitude  and  phase  of 
SAWs  in  the  product  being  tested. 

Measurement  of  the  phase  of  SAWs  makes  significant  demands  on  the  stability  of 
the  electrical  lengths  of  the  measurement  and  reference  channels  of  the  installation. 
The  requirements  are  considered  by  equating  the  optical  lengths  of  the  channels 
after  the  acoustico-optical  modulator,  and  through  construction  measures  which 
stabilize  and  minimize  them.  At  the  same  time,  adjustment  mirror  7  provides 
colinearity  of  the  reference  ray  in  the  forward  and  reverse  directions. 

Due  to  the  smallness  of  the  diffraction  angles  in  the  installation,  the  divi.sion  of  the 

orders  of  diffraction  is  done  by  a  prism  with  two 
reflecting  surfaces  4.  This  construction  significantly  decreases  the  effect  of 
adjustment  error  in  the  reflecting  elements  on  the  accuracy  of  bringing  the  rays  to 
the  receiving  window  of  the  photodetector. 

The  described  structure  of  the  optical  sensor  makes  it  possible  to  .study  SAWs  in  a 
wide  range  of  frequencies. 

A  .significant  element  of  the  optical  sy.stem  of  the  installation  is  the  acoustico- 
optical  modulator.  At  low  frequencies  a  liquid  Bragg  cell  with  an  excitation 
frequency  of  28±5  MHz  is  used  for  this  modulator.  In  terms  of  construction,  the 
cell  is  in  the  form  of  a  cuvette  made  of  a  material  with  significant  acoustic 
ab.sorption.  The  cuvette  is  filled  with  distilled  water.  A  removable  piezoelectric 
transformer  is  attached  to  the  cuvette  in  the  window  with  a  water-tight  seal. 
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A  mode  of  traveling  acoustic  waves  is  provided  by  choosing  a  sufficiently  long  cell 
and  by  using  a  final  acoustic  absorber.  The  size  of  the  cell  in  the  direction  of  the 
optical  beam  (L)  provides  operation  in  a  Bragg  diffraction  mode: 


where  A  is  the  length  of  the  acoustic  wave  in  the  liquid  filling  the  cell  and  is 
the  wavelength  of  the  laser  radiation. 

In  the  cell  used,  L  =  60  mm,  which  provides  a  Bragg  diffraction  mode  at  frequencies 
greater  than  10  MHz.  At  frequencies  greater  than  60  MHz  solid  state  acoustico- 
optical  modulators  are  used. 

For  fine  tuning,  the  cell  is  placed  at  an  adjustment  node.  The  criterion  of  optimal 
adjustment  is  the  maximum  of  the  signal  at  the  output  of  the  photodetector  (the 
adjustment  limit  is  ±0.5**).  There  are  analogous  adjustment  nodes  at  the  mirror  of 
the  reference  channel  and  at  the  object  stage  10  on  which  the  product  to  be  studied 
is  placed.  The  linear  mechanical  sweep  of  the  stage  is  transformed  into  an 
electrical  signal  which  enters  input  X  of  the  PDS-021  automatic  recorder. 

Information  on  the  amplitude  and  phase  of  SAWs  is  separated  out  by  processing  the 
signals  of  the  photodetector  with  an  electronic  unit,  which  is  a  two-channel 
superheterodyne  receiver  with  a  cascade  of  automatic  amplification  regulation  which 
is  common  to  both  channels.  Two  signals  are  processed  in  the  receiving  unit  (1): 
2/s  and  2/b  +  /c.  A  feature  of  the  schematic  is  the  use  of  an  additional  mixer  in 
the  second  channel  to  decrease  the  frequency  of  the  (2/^  +  /c)  signal  to  the 
frequency  of  the  2/b  signal,  as  well  as  the  use  of  the  same  type  of  intermediate 
frequency  amplifiers  at  2/^  in  both  channels. 

The  installation  uses  acoustico-optical  modulators  at  the  following  frequencies: 
28 ±5  MHz,  80 ±1  MHz,  and  120J::3  MHz.  Depending  on  the  frequency  of  the 
modulator  the  necessary  intermediate  frequency  amplifiers  are  installed  in  the 
receiving  electronic  unit.  The  output  signals  of  the  channels  are  reduced  to 

465  kHz  by  reducing  mixers  and  the  appropriate  heterodynes  in  the  intermediate 
frequency  amplifier  complex.  These  signals  enter  the  input  of  a  standard 

phasometer,  which  measures  the  phase  of  the  SAWs.  The  signal  of  the  amplitude  of 
the  SAWs  is  the  change  in  the  control  voltage  of  the  automatic  amplification 
regulation  cascade  at  the  input  of  the  receiving  electronics  unit.  The  amplitude  and 
phase  profiles  of  the  SAWs  are  recorded  in  order  by  an  automatic  recorder. 

Figures  2  and  3  show  the  amplitude  and  phase  profiles  of  SAWs  measured  in  a  filter 
with  fo  =  28  MHz.  From  the  amplitude  profile  (Figure  2a)  and  the  noise  curve 

(Figure  2b)  of  the  measurement  channel  one  can  conclude  that  the  upper  curve 

transmits  the  deformation  of  the  surface  by  SAWs  of  the  studied  filter  with  a  high 
degree  of  accuracy.  It  is  obvious  that  surface  oscillations  As  are  grouped  toward 
the  longitudinal  axis  of  the  product.  The  phase  profile  (see  Figure  3)  makes  it 
possible  using  known  means,  to  determine  the  rate  of  propagation  of  the  SAWs,  The 
described  installation  is  intended  for  the  study  of  SAW  parameters  in  the  range  of 
frequencies  from  15  MHz  to  300  MHz. 

The  authors  thank  A.  D.  Belyayev  for  a  number  of  valuable  comments  expressed 
during  discussion  of  the  work. 
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Figure  2  [left].  Amplitude  distribution  of  SAWs  in  the  transverse  cross  section  of 
the  sample  (a)  and  the  noise  curve  (b)  of  the  measurement  channel,  c.  As,  relative 
units;  d.  aperture  of  the  transformer. 


Figure  3  [right].  Distribution  of  SAW  phase  in  the  direction  of  propagation.  X-axis 
in  Aim, 
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[Text]  Providing  electromagnetic  compatibility  to  electrical  power  consumers  with 
equipment  having  nonlinear  volt-ampere  characteristics  and  supplying  them  from  the 
power  grid  requires  the  solution  of  a  number  of  complex  problems,  including  that  of 
calculating  non-sinusoidal  modes  for  the  electrical  supply  system  The  theoretical 

and  experimental  studies  we  have  carried  out  on  operating  modes  for  higher-harmonic 
sources  have  shown  that  over  long  periods  of  time  the  higher-harmonic  currents  they 
generate  vary  in  a  random  way  due  to  the  stochastic  nature  of  variation  in  the  load 
over  these  devices  and  the  voltage  of  the  power  grid.  The  ever-changing 
configuration  of  an  industrial  power  grid  and  changes  in  the  parameters  of  its  elements 
are  also  random  in  nature.  In  this  regard,  the  development  of  a  probabilistic- 
statistical  method  for  calculating  higher-harmonic  modes  has  significant  practical 
value.  In  this  article  we  describe  a  probabilistic  method  we  have  developed  for 
calculating  the  higher-harmonic  modes  of  an  industrial  power  supply  system  when  the 
mode  parameters  of  its  sources  vary  in  a  stochastic  manner. 

Calculating  a  higher-harmonic  mode  includes  determining  the  amplitude  spectrum  of  the 
node  voltages  and  currents  in  the  loops  of  the  supply  system  equivalent  circuit,  and 
the  coefficient  of  higher-harmonic  content  (non-sinusoidality)  of  the  voltage  and 
equivalent  effective  higher-harmonic  currents  in  the  electrical  grid  elements.  The 
stochastic  nature  of  variation  in  the  parameters  of  a  higher-harmonic  source  leads  to 
random  variation  in  these  higher-harmonic  mode  parameters.  In  this  regard,  the 
mathematical  model  which  most  adequately  describes  a  non-sinusoidal  mode  is  a 
probabilistic  model.  The  parameters  of  the  grid  elements  and  their  configurations  are 
considered  constant  and  are  introduced  into  the  calulation  via  a  determinant  operator, 
while  the  nonlinear  load-induced  higher  harmonics  are  modelled  with  random  variables 
having  correlated  real  and  imaginary  parts.  For  the  majority  of  practical  problems  it 
is  wholly  satisfactory  to  limit  oneself  to  a  determination  of  the  numerical  probabilistic 
characteristics  of  the  higher-harmonic  mode  parameters  within  the  framework  of  a 
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correlational  theory — one  with  mathematical  expectations,  dispersions  and  correlational 
moments.  A  method  for  numerical  calculation  of  higher-harmonic  mode  parameter 
characteristics  in  terms  of  the  numerical  characteristics  of  the  higher-harmonic 
current  source  is  presented  below. 

A  higher-harmonic  mode  calculation  consists  of  the  following  stages: 

1.  Construct  an  equivalent  circuit  for  the  higher-harmonic  power  supply  system. 

2.  Calculate  the  equivalent  circuit  element  resistances  at  the  higher-harmonic 
frequencies  in  question. 

3.  Determine  the  probabilistic  characteristics  of  the  higher-harmonic  current 
generated  by  devices  having  nonlinear  volt-ampere  characteristic  curves. 

4.  Consolidate  all  higher-harmonic  sources  connected  to  the  same  node  of  the 
electrical  circuit. 

5.  Calculate  the  numerical  characteristics  of  the  node  voltages  and  higher-harmonic 
currents  in  the  loops  of  the  power  supply  system  equivalent  circuit. 

6.  Determine  the  numerical  characteristics  of  the  coefficient  of  non-sinusoidality  of 
the  node  voltages  and  equivalent  effective  values  of  the  higher-harmonic  current  in 
the  loops  of  the  equivalent  circuit. 

The  first  two  stages  are  performed  by  a  common  method  used  for  deterministic 
calculations  of  higher-harmonic  modes,  so  we  will  not  look  into  that  in  this  article. 

The  third  stage  includes  a  calculation,  for  each  harmonic,  of  the  mathematical 
expectations  of  the  real  i/])  and  imaginary  i/])  components  of  the  random 

current  vector  Jp,  the  dispersion  of  Jp  along  the  coordinate  axes  DlJ'p]  and  DU  p], 
and  the  correlation  coefficient  of  the  real  and  imaginary  components  of  Jp.  For 
the  canonical  harmonics  of  rectifier  converters  and  thyristor  voltage  regulators 
commonly  used  in  industry,  these  parameters  may  be  calculated  by  the  following 
formulas: 


1.  Controllable  and  noncontrollable  rectifiers  with  a  small  range  of  variation  in  the 
angle  of  regulation  a.  (rectifiers  for  shop  dc  circuits,  converters  for  ganged  dc  electric 
motors,  non-reversible  converters  for  electric  motors  with  velocity,  voltage  or  motor 
emf  feedback,  etc.): 

mJ'p]  -=AdUp]cos{lMm>, 

AllJp]  =A4lJp]sini'^A4m); 

D{J'A  =DU  i']cos^(7A/[^J); 

DU'p]  =D[J..]sin"(7MV)]). 

For  7: 

=  if  MU' pUIU"  p\>0-, 

^bmp~  0.5,  if  MU  p]MU  ‘p]<^0, 

For  13: 


BM./' 


=0. 


2.  For  converters  using  current  stabilizers  and  three-phase  thyristor  voltage 
regulators  for  metallurgical  electric  furnaces: 

MU'p]  =A/IJ..](l  -^DliD])cos(i/Ml£i]); 


MU'p]—MUp]^1  — ^Df?l>]]sin(7Af[il)]); 


14 


D\J'A=u^M^lJADmsinHuA4my, 

_  _  Al^lJADU:-]  f/M^[  JllDlJl'l 
®"'"  IMUMIUM  DU'MJu] 


3.  Converters  with  a  large  range  in  variation  of  the  angle  of  regulation,  reversible 
and  non-reversible  electric  motor  converters,  converters  for  plasma  steel  smelting 
furnaces,  vacuum  furnaces,  etc.)  and  three-phase  thyristor  voltage  regulators  for 
deposition-type  electric  furnaces: 

/W[<7 

D[J' J  =D[  V 

The  parameters  AdlJv]  and  D[JA  entering  into  the  formulas  are  the  mathematical 
expectation  and  dispersion  of  the  modulus  of  the  harmonic  of  a  converter,  and 
MW  and  DW  are  the  mathematical  expectation  and  dispersion  of  the  angle  +7/2 

are  defined  by  the  following  expressions: 


1)  for  converters  with  a  high  range  of  variation  in  the  angle  of  regulation, 


where 


Here  is  the  nominal  rectified  current  of  the  converter; 

2)  for  converters  with  a  small  range  of  variation  in  the  angle  of  regulation. 


2fei^„Aiirj 

xi/lW[7] 


|sin(^yW[7ni; 


D[J.]==^DU*J; 

/WI^]=Adla}+0.5./W[7]; 


where 


DW  -  +DlU%]cot"M[oL], 

A/[al  =arccos{Af[(./dl  +0.5xtM[Id\h 


Af  [7] =arccos{cos  Af(a] — x*M[l*d])  —  Af[a]. 


Here  t/J  and  Ul  are  the  relative  rectified  voltage  and  the  grid  voltage,  is  the 
relative  inductive  reactance  in  the  commutator  circuit  of  the  transformer,  and  a  and 
7  are  the  angles  of  regulation  and  commutation; 

3)  for  noncontrollable  rectifiers 


AAIJJ 


ir//(i/-l)M7ol 


|sin(^=^-^A/[7oJ)l; 


+ft0l=|Afl7ol. 
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where 


M7o]=arccos{1  -xtM[I%]) 


4)  for  three-phase  thyristor  voltage  regulators 


Dti&I  = 


Mm- 

4.4 

M-[Uc^ 


.7-2.1 


(Wt/, 


M[Uc]  ’ 


out 


1 


MWc] 


Here  is  the  regulator  load  current  for  a=0,  and  is  the  effective  output 
voltage  of  the  regulator. 


Consolidating  equivalent  higher-harmonic  sources  is  done  when  several  higher- 
harmonic  sources  are  connected  to  the  same  node  of  an  electronic  circuit. 

The  equivalent  source  should  also  have  a  probability  characteristic  such  that,  in 
calculating  the  harmonic  mode  parameters,  the  source  mode  parameters  will  be  the 
same  before  and  after  consolidating  the  sources. 


A  theoretical  investigation  has  shown  that  the  distribution  of  the  equivalent  random 
vector  Jut:  becomes  asymptotically  closer  and  closer  to  a  bivariate  normal  distribution 
as  the  number  of  harmonic  sources  joined  together  increases.  Its  numerical 
characteristics  may  be  determined  from  formulas  obtained  by  using  the  rule  for 
summing  independent  random  vectors: 


m 


k  =  l 


^BMa/ 


2M[j'uWy'M  Dij'ADm 


For  minimizing  the  higher-harmonics  content,  we  need  to  know  the  probabilistic 
characteristics  of  the  modulus  of  J Theoretical  and  experimental  studies  we  have 
performed  showed  that  the  distribution  of  the  modulus  of  the  equivalent  current  J 
is  close  to  a  normal  distribution,  but  with  a  positive  asymmetry  that  decreases  as  the 
correlation  coefficient  of  the  components  of  J increases.  Decrease  in  the 
coefficient  of  variation  in  the  J components  also  brings  the  distribution  equation 
for  the  modulus  of  J closer  to  normal.  So  substituting  in  a  normal  distribution  for 
the  the  current  modulus  distribution,  we  obtain  the  following  expression  for 
determining  the  mathematical  expectation  M[J^y]  and  dispersion  DU „y]  of  the  modulus 

cl*  fZ:'  4  p- _ 

^|Af^l42:]-0.5D[J^r];  (1) 


16 


(2) 


where 


D\<J j/i;] =A/[t7 27^/1  —  [«7 j/^]  0.5D[t7 t/z;]f 

Aj\  jf. .]  =  Af  JU]  +AriJ:r]-\-D[J,r:]  +D\J:.r\; 

DUl^]  =4.4^ J;,r]  +4M1  J';r]D[  Jlirl +2D-V„7:1 +2D"[  + 


Formulas  (1)  and  (2)  are  obtained  for  the  general  case  in  which  the  components  are 
correlated,  dispersion  along  the  axes  is  not  equal,  and  the  mathematical  expectation  is 
non-zero.  In  the  particular  case  of  a  node  being  connected  to  a  higher-harmonic 
source  with  uniformly  distributed  phase  in  the  range  from  0  to  27r,  the  numerical 
characteristics  of  the  equivalent  current  modulus  may  be  defined  in  a  very  simple 


equation: 


■It 


//c 


k=l 


(3) 


which  is  valid  for  connected  sources.  Here  the  current  modulus  distribution  is 

subject  to  the  Rayleigh  equation.  The  maximum  error  in  formulas  (3)  and  (4),  in 
comparison  to  (1)  and  (2),  is  no  more  than  13.5%. 

!>,,,,  the  voltage  of  the  harmonic  at  the  power  supply  system  equivalent  circuit 
node,  is  a  linear  combination  of  the  driving  currents  of  the  nodes 

77? 

k^l 

Here  is  the  complex  nodal  resistance  between  the  and  fc*'"  nodes,  is  tl^e 

complex  random  current  from  the  harmonic  of  the  source  connected  to  the  k 
circuit  node,  and  m  is  the  number  of  nodes  connected  to  the  harmonic  source. 

Ui,i  is  completely  described  by  the  simultaneous  bivariate  distributions  of  a  system  of 
random  quantities  describing  its  real  U'^i  and  imaginary  U'U  components.  Here  then 

MVL},^  =  M[UU]  +jM[U:,]; 

DlO,i]  =  DUr,t]+jDlLni]. 

Using  expression  (5)  and  the  theorem  about  the  numerical  characteristics,  we  obtain 
the  following  expression  for  the  mathematical  expectation  and  dispersion  of  the 
random  vector  components: 


M [U'v,]  =  H  [Tv.]  +  X,  [T'v/:]  ; 

-  1  ^  I 

M [ U"vi]  =  Ij  ’ 


(6) 

(7) 
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(8) 


DiU\i]  =  + 

ft-l  /t-1 

m 

+  2  Rvih^  vih}^D  £)  ^ BMvft  J 

/i=l 

nj.  ?n 

WZ)[/'v/J+X  /?w/r^[/"v;.]  + 

fc-i  ft=i 

m 

+  2  R^i„X,,niD[J\u]D[J\u]  r„,vft.  (9) 

/l=I 

Here  is  the  correlation  coefficient  of  the  real  and  imaginary  components  of 

For  solution  of  practical  problems  involving  determining  the  coefficient  of  non- 
sinusoidality  of  the  voltage,  calculation  of  damage  caused  by  higher  harmonics,  and 
choice  of  filter  parameters,  we  need  to  know  the  numerical  characteristics  of  the 

random  modulus  of  the  voltage  of  a  harmonic  at  the  node  ,  which  is 

completely  defined  by  a  univariate  distribution.  For  determining  the  mathematical 

dispersion  of  the  modulus  of  the  voltage  of  the  harmonic  at  the 
i  "  circuit  node,  we  obtain  the  following  expression: 

M [ f/,,]  =  yAP [ ]+M^[ U%i]  +  D  [ U\,]  +  D  [ U%,- ]-D[ (10) 

D[Uvi]  ^D[U^i]  —  A^;  (11) 

m 

K  [  U\.i  U'\.,]  =  ij],  ^v,h2fvWi  {D  [/'vft  ]  -  D  [/"v„ ] )  + 

/i-I 

m 

+  H  reMvA.  (  )  2) 

The  right-hand  side  of  equation  (11)  contains  two  components:  b[U^i],  the  dispersion 
component  obtained  by  ^linearizing  the  function  U^i  in  the  neighborhood  of  the  point 
and  A^,  a  correction  for  the  nonlinearity  of  L/,.;.  By  virtue  of  the 

expression 


5[(/v,]  = 

A12  [  f/'v,]  D  [  U\i]  +  [  U'\i]  D  [  U\i]  +  2/Vf  [  U\i]  M  [  U\i]K[U\i  (/%,] 


(13) 


A2  = 


{/W2 [(/"v,] D  [ U\i]  +  [ U\i] D [ U\i\  -2M[U\i]  M [ U\i]  K [ (/%■] }2 


A{iW[U',i]+iW[U\,'\Y 


(14) 


The  numerical  experiment  we  conducted  confirmed  the  need  for  the  nonlinearity 
correction  A  .  E.g.,  for  T'gi.,  =0  and  the  coefficients  of  variation  in  the  components  of 
the  vector  U equal  to  two,  without  the  correction  A^  the  error  in  determining  AIIU.,] 
is  -1-15%,  and  for  is  —50%. 
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The  harmonics  current  in  the  loops  of  the  equivalent  circuit  is  a  linear 
combination  of  the  driving  currents  in  the  nodes: 

m 

T.  J  i/k  • 

Here  C'^r,,=-C'^jk+jC'f,ii.  is  the  complex  coefficient  of  the  current  distribution  for 
the  loop. 

The  equations  (5)  and  (15)  are  identical  in  structure.  Therefore  the  expressions  for 
determining  the  mathematical  expectation  and  dispersion  of  the  components  of  the 
vector  will  be  similar  to  expressions  (6)  through  (9),  but  instead  of  real  and 
imaginary  components  of  nodal  resistances,  we  need  to  insert  the  real  and  imaginary 
parts  of  the  coefficients  of  the  driving  current  distribution  C^jk.  The  mathematical 
expectation  and  dispersion  of  the  modulus  of  1^,;  may  be  calculated  from  the  formulas 
(10)  through  (14),  in  which  instead  of  the  mathematical  expectation  and  dispersion  of 
the  components  of  we  need  to  insert  the  mathematical  expectation  and  dispersion 
of  the  components  of  the  random  vector  Ij,j. 

The  coefficient  of  non-sinusoidality  of  the  voltage  in  the  nodes  of  the  supply 
system  is  a  nonlinear  function  of  the  random  moduli  of  the  higher-harmonic  voltages 
and  the  voltage  at  the  fundamental  frequency.  To  a  first  linear  approximation,  the 
mathematical  expectation  and  dispersion  may  be  calculated  by  the 

formula 


MKhJ-  If 


D[Xhc1 


irr 

ul 


n 


(16) 


(17) 


Here  n  is  the  number  of  higher  harmonics,  and  (/„  is  the  nominal  grid  voltage. 

Using  the  formulas  (16)  and  (17)  is  permissible  for  coefficients  of  variation  in  the 
voltage  harmonics  of  no  more  than  0.5  and  positive  correlation  coefficients  between 
the  voltage  harmonics  in  the  nodes. 

The  error  in  determining  here  is  no  more  than  10%.  For  larger  coefficients  of 

variation  we  need  to  include  a  correction  for  the  nonlinearity  in  the  function 

The  numerical  characteristics  of  the  equivalent  effective  current  in  the  supply 
system  elements  are  determined  from  the  equations 


Mile] 


D[/e]=  - 
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The  underlying  theoretical  hypothesis  of  this  calculation  method  is  confirmed  by  the 
results  of  experimental  investigations  conducted  in  ceramics  and  petrochemical 
industrial  complexes.  The  deviation  between  the  calculated  and  experimental  higher- 
harmonic  mode  parameters  is  not  more  than  ±30%. 

Conclusions.  1.  In  this  article  we  have  proposed,  in  solving  the  question  of 
electromagnetic  compatibility,  to  describe  non-sinusoidal  mode  parameters  using  random 
models  and  systems  of  random  quantities  with  correlated  real  and  imaginary 
components.  In  the  majority  of  cases,  determining  just  the  numerical  probabilistic 
characteristics  of  higher-harmonic  mode  parameters  is  enough  for  solution  of  most 
practical  problems. 

2.  A  probabilistic-statistical  method  has  been  developed  for  calculating  the  higher- 
harmonic  modes  which  enables  us  to  determine  the  mathematical  expectations  and 
dispersions  of  higher-harmonic  voltages  in  the  nodes  and  harmonic  currents  in  the 
loops  of  the  supply  system  equivalent  circuit,  and  the  probabilistic  numerical 
characteristics  of  non-sinusoidality  coefficient  of  the  voltages  and  the  equivalent 
effective  currents. 
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[Article  by  V.A.  Rodionova  and  N.A.  Fomina] 

[Text]  Plate  and  shell-type  structural  elements  made  of  polymer  composite  materials 
have  become  commonly  employed  in  designs  over  the  past  few  years,  and  this  has  led 
to  the  need  for  a  theory  of  anisotropic  plates  and  shells  that  would  take  into  account 
the  specific  characteristics  of  composite  materials,  such  as  their  low  compliance  to 
shear,  and  low  strain  resistance  in  the  direction  normal  to  the  coordinate  surface. 

A  (modified)  theory  of  anisotropic  plates  was  constructed  in  *,  which  may  be  viewed  as 
a  first  approximation  when  solving  a  boundary  problem  in  a  three-dimensional  theory 
of  elasticity  using  the  weighted  residuals  method.  When  employing  this  method,  static 
boundary  conditions  over  the  faces  of  the  plate  and  the  Saint  Venants  condition  over 
the  cylindrical  side  surface  of  the  plate  are  exactly  satisfied.  It  would  be  of  some 
interest  to  construct  a  second  approximation  (a  refined  modified  theory)  that  would  in 
addition  take  into  account  certain  surface  load  components,  self-compensating  over  the 
thickness  of  the  plate,  that  act  on  the  side  surface  of  the  plate.  It  is  the 
construction  of  just  such  a  theory  that  is  to  be  examined  in  this  article. 

The  proposed  theory  may  be  constructed  by  hypothesis.  Below,  we  will  use  this  very 
approach  for  construction  of  a  theory  of  plates  corresponding  to  a  second 
approximation  using  the  weighted  residuals  method. 

In  order  to  make  the  situation  clear,  we  will  limit  ourselves  to  looking  at  constant¬ 
thickness  anisotropic  plates  in  the  absence  of  body  forces.  A  plate  will  be  viewed  as 
a  thin  anisotropic  elastic  solid  possessing  a  plane  of  elastic  symmetry.  An  orthogonal 
curvilinear  coordinate  system  ai,  and  z  is  associated  with  the  midplane  of  the 
plate,  with  the  linear  coordinate  z  directed  perpendicular  to  it.  The  primary  direction 
of  elasticity  coincides  with  this  linear  coordinate. 

Static  boundary  conditions  are  assigned  over  the  face  plane  of  the  plate  (z  =  ±h/2) 

o^3=qrt,  0^3=02.  o^3=<7^  for  T  =  ±h/2  (1) 

(where  h  is  the  plate  thickness).  We  introduce  the  following  hypotheses: 

— (statistical):  The  perpendicular  tangential  strains  a  ,3  and  0^23  and  normal  strain  0-33 
are  distributed  over  the  plate  thickness  according  to  polynomial  equations  that  are, 
respectively,  of  fourth  and  fifth  order  in  z,  i.e* 


(kinematic):  the  tangential  Ui  and  xt2  and  normal  ^^3  components  of  the  displacement 
vector  are  distributed  over  the  plate  thickness  according  to  polynomial  equations  that 
are,  respectively,  of  fifth  and  fourth  order  in  z,  i.e. 


A=0 


.9  =  0 


or 


«.=:  C),  y 


(3) 


In  these  assumed  expressions,  crj3,(733,  and  are  functions  of  the  coordinates  and 
0^21  are  Legendre  polynomials,  qi=qt  —QV,  and  m;  =(9^ 4-gf )?i/2. 

Integrating  expression  (2)  with  weighting  by  P,:„  Pi  and  P2  with  respect  to  z  from 
z  =  —h/2  to  z—  -\-h/2,  we  find  that 
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Introducing  the  integral  quantities 
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and  using  the  third  equilibrium  equation  (21.5)  from  ^  into  which  is  substituted 
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Henceforth  we  will  use  the  following  representations  for  the  stresses  and  strains: 

a,..=  Vo|))P,(C;,  .,.=  ;Vg<5)P,(C);  /.  y  =  l,  2; 
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/2=J 
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^,,3  =  V  [D ,  =  ^  e^P,  (Q. 
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By  directly  integrating  with  respect  to  z  from  z  =  — h/2  to  x—h/2  with  weighting 
factors  P,-,  using  equations  (21.1)  from  ^  we  find  the  dependence  between  the  strain 
components  cTj'  and  the  displacements  u,  .  In  a  similar  manner,  we  find  the  Hookes 
Law  relations  using  expressions  (2.5.18)  and  (2.5.19)  from  . 


In  equation  (3)  there  are  seventeen  unknown  functions  u,,  representing  all  the 
components  of  the  stress-strain  state  of  the  plate.  By  using  equation  (2),  obtained  by 
consideration  of  boundary  conditions  (1)  over  the  face  surface  of  the  plate,  we 
express  unknown  functions  (i=l,  2;  fc=4,  5)  and  W3  (fc=3,  4)  in  terms  of  the 
others.  Then,  instead  of  seventeen,  there  are  only  eleven  unknown  Ui  . 


We  introduce,  in  accordance  with  the  introduced  hypothesis,  the  following  expression 
for  the  potential  energy  of  plate  deformations 


v=^ Jj  s +  T + T 


5  \  i, 

9 


A-  -L  -j-  2-  "V 

i  7  i;  u )  '  9 


We  obtain  equilibrium  equations  for  the  plate  elements  and  the  natural  boundary 
conditions  using  the  Lagrange  variational  equation  SV  —6 (in  the  absence  of  body 
forces),  in  which  the  work  of  external  surface  forces  over  the  corresponding  virtual 
displacements  6A^6A+^8A-+SA,  is  equal  to  the  work  of  surface  forces  applied  to 
the  upper  (z=h/2)  and  lower  (z  =  — h/2)  face  planes  of  the  plate,  as  well  as  to  the 
cylindrical  side  surface  S*. 

One  may  show^  that  the  system  of  partial  differential  equations  we  are  solving  is 
twenty-second  order  overall.  In  this  case  that  system  is  broken  down  into  two 
systems  of  equations  of  tenth  and  twelfth  order,  which  describes  respectively 
symmetric  and  antisymmetric  stress-strain  states  of  the  plate.  On  each  edge  of  the^ 
plate  there  are  eleven  boundary  conditions.  This  question  is  examined  in  detail  in  . 

Comparison  of  the  results  of  calculations  with  this  theory  for  an  isotropic  cantilever 
beam  with  transverse  tangential  load  applied  to  the  free  end  in  a  quadratic-law 
distribution  over  the  beam  thickness^  with  the  results  of  solution  of  the  apjropnate 
boundary  problem  in  the  plane  elasticity  theory  by  the  method  proposed  in 
indicates  that  the  proposed  theory  and  its  basis  enables  us  to  find  the  secondary 
Saint-Venant  edge  effect-type  stress-strain  states  corresponding  to  the  first  of  the 
two  roots  of  the  transcendental  equations  (r)  and  (w)  in  '  ”  ,  i.e.  it  allows  us  to  take 

into  account  the  effect  of  these  surface  load  components,  self-compensating  over  the 
plate  thickness  and  applied  to  the  side  surface  of  the  plate,  which  vary  according  to  a 
polynomial  of  no  higher  than  third  order  in  z. 
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[Article  by  Yu.  V.  Yadykin,  V.  A.  Gorban,  and  V.  M.  Kolotovich,  Kiev] 

[Text]  /The  dynamics  of  low-frequency  vibration  protection  systems  which 
contain  a  compensator  of  imposed  force,  along  with  elastic  and  damping 
elements,  is  investigated.  Results  are  given  from  numerical  analysis  of 
several  structural  layouts  of  vibration  protection  systems.  It  is  shown  that 
vibration  protection  systems  built  on  the  proposed  principle  possess  a  broader 
working  frequency  range  of  vibration  protection./ 

With  the  development  of  ocean  exploration  facilities,  the  issue  of  vibration 
protection  of  highly  sensitive  apparatus  against  the  action  of  low-frequency 
vibrational  loads  generated  by  the  rolling  motion  of  the  water  is  becoming 
increasingly  urgent.  Traditional  methods  of  vibration  protection  involving 
the  use  of  passive  elements  (inertial,  elastic  and  dissipative)  in  many  cases 
prove  ineffective  [1-3]. 

Analysis  of  the  amplitude -frequency  characteristics  of  vibration  protection 
systems  [3-5]  suggests  one  major  conclusion- - in  order  to  expand  the  working 
frequency  range  it  is  necessary  to  reduce  the  resonance  frequency  of  the 
system,  which  exerts  major  influence  on  the  lower  limit  of  the  working  range. 
This  can  be  achieved  by  lessening  the  rigidity  of  the  elastic  element  or 
increasing  the  mass  of  the  body  being  stabilized,  which  is  not  always 
practicable.  The  change  in  rigidity  of  the  elastic  element  and  increase  in 
mass  of  the  stabilized  body  are  limited  by  the  permissible  static  deformation 
of  the  elastic  element  under  the  action  of  the  stationary  load  on  the  body. 

For  a  limited  stroke  of  the  elastic  element  5  <  A  •=  F/C,  the  frequency  of 
natural  oscillations  of  the  stabilized  body  will  be: 


/7,.=1  Cpn^y  Flm!^ 

and  depend  on  the  mass  of  the  body  m  and  the  permissible  deformation  of  the 
elastic  element. 

It  is  not  difficult  to  see  that  the  requirements  imposed  on  vibration  protec¬ 
tion  systems  are  contradictory  [6] .  On  the  one  hand,  the  elastic  element  of 
such  systems  should  be  compliant,  in  order  to  achieve  a  low  resonance 
frequency;  on  the  other  hand,  it  should  be  quite  rigid,  in  accordance  with  the 
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conditions  of  static  equilibrium.  This  narrows  the  range  of  effective  use  of 
vibration  protection  systems. 

In  the  present  article,  we  set  forth  a  principle  for  design  of  autonomous 
systems  of  low-frequency  vibration  protection- -low-frequency  decoupling  of 
both  passive  (no  energy  supplied)  and  active  type,  based  on  introducing  a 
stationary  force  input  in  the  compensator  system.  We  also  present  results  of 
an  investigation  of  the  amplitude -frequency  characteristics  of  low-frequency 
decoupling  for  ocean  exploration  facilities.  Introduction  of  the  compensator 
makes  it  possible  to  eliminate  the  basic  engineering  conflict  of  vibration 
protection  systems- -significantly  reducing  the  stiffness  of  the  elastic 
element  of  the  vibration  insulator  without  changing  its  static  range.  It  is 
shown  that  effective  low-frequency  decoupling  has  an  elastic  element  of  low 
rigidity.  In  such  case,  the  length  of  coupling  of  the  source  of  perturbation 
and  the  body  being  stabilized  is  self -regulating,  i.e.,  change  occurs  without 
employing  automatic  control  systems  as  a  function  of  the  relative  displace¬ 
ments  of  the  source  and  the  body  being  stabilized. 

In  the  coupling  that  joins  the  body  to  the  source  of  perturbation,  let  us 
divide  the  force  F  into  static  and  dynamic  F^j  components.  The  static 
component  Fg^  will  be  compensated  by  a  force  of  constant  magnitude  and 
direction,  independent  of  the  velocity,  the  acceleration,  and  the  relative 
position  of  the  body  being  stabilized,  while  the  influence  of  the  dynamic 
component  F^  on  the  movement  of  the  body  being  stabilized  will  be  diminished 
by  altering  the  length  of  the  coupling. 

Among  the  familiar  forces  in  mechanics,  we  may  identify  those  which  satisfy 
the  requirement  of  constant  magnitude  and  direction  for  the  feasible  movements 
of  the  object  and  the  source  of  the  perturbations.  Such  forces  include,  for 
example,  the  force  of  gravity,  the  Archimedean  force,  the  centrifugal, 
hydrodynamic  force  of  a  screw,  and  so  on. 

The  layout  for  a  passive  low-frequency  vibration  protection  system  includes 
the  object  (body  being  stabilized)  and  the  source  of  perturbations,  joined 
together  by  a  low-frequency  decoupling,  in  which  a  compensator  of  the  imposed 
force  is  placed  in  parallel  with  a  vibration  insulator.  The  amplitude - 
frequency  characteristic  of  the  stabilized  body,  joined  to  the  source  of 
perturbations,  has  the  same  appearance  as  an  object  connected  to  a  source  of 
perturbations  by  an  elastic  element  in  the  trivial  layout  [2,  6].  However, 
the  resonance  frequency  for  the  system  (Fig.  1)  is  shifted  into  the  region  of 
low  frequencies  of  oscillation,  i.e.,  po*  =  J^/m  =  jF^/mA  and  tends  toward 
zero  at  the  limit. 

As  an  example  of  the  implementation  of  the  principle  of  construction  of  a  low- 
frequency  decoupling,  let  us  examine  several  structural  layouts  for  vibration 
protection  systems  of  this  type. 


Figure  1  shows  a  layout  in  which  the  force  of  buoyancy  or  gravity  is  used  as 
the  compensating  static  load.  A  body  of  mass  m^  is  joined  to  the  source  of 
perturbations  across  a  rocking  lever  with  ideal  hinge  by  means  of  a  flexible 
nonstretching  string,  passed  across  the  hinge.  At  the  free  end  of  the  string 
is  suspended  an  auxiliary  body  (compensator)  of  mass  m2.  Both  bodies  are 
joined  to  the  rocking  lever  by  elastic  elements.  The  system  is  in  static 
equilibrium.  The  hinge  enables  self  regulation  of  the  length  of  the  coupling 
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between  the  body  and  the  source  of  perturbations  (the  length  of  the  segment  of 
the  coupling  between  the  auxiliary  body  and  the  source) .  In  the  case  when  the 
stabilized  body  and  the  auxiliary  body  are  placed  in  a  liquid,  the  equation  of 
motion  of  the  system  is : 

m:x,+R,\x,\x-C,{xo-x,)=m.;x^+R^\x■M^-C^{x-x^),  (1) 

where  x^  +  X2  =  2xo,  Xq  -  h  sin  flt  is  the  external  perturbation;  =  O-SC^^ipSi, 
R2  =  0.5Cx2PS2;  mj^*,  m2*  are  the  masses  of  the  stabilized  body  and  the 
compensator  (auxiliary  body) ,  taking  into  account  the  added  masses  of  the 
liquid;  R^,  R2  are  the  parameters  of  hydrodynamic  resistance  of  the  body  and 
the  compensator;  C^,  C2  are  the  rigidities  of  the  elastic  elements;  C^2 

are  the  coefficients  of  hydrodynamic  resistance;  Sj^,  S2  are  the  cross  section 
areas  of  the  stabilized  and  auxiliary  bodies  at  the  maximum  midsection;  p  is 
the  density  of  the  medium. 


Fig.  1.1-  body  being  stabilized,  2  - 
source  of  perturbations ,  a  rocking  lever 
with  ideal  hinge,  3  -  flexible  nonstretch¬ 
ing  string,  4  -  auxiliary  body,  5,6- 
elastic  elements. 


Using  the  nondimens ional  variables: 


(2) 


we  convert  equation  (1)  to  the  form: 


where 


|+/?/V|g|i/(M+l)+iV12  cos  T-sl§/ {M+i)+llq-= 
=  2A^|2  cos  T— Ijcos  ]sin  t, 
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(3) 
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(4) 


Using  the  above  definitions,  we  find  for  the  tension  T  of  the  string: 
o  - ^ 1  +  [m:Qi+H,hQ^  I  §  I  i+C,  (5-sin  x)  ] , 

mt  g  Till  g 

where  nij^'  =  is  the  mass  of  liquid  displaced  by  the  body  being 

stabilized. 

Results  of  calculations  of  the  d3mamic  characteristics  (amplitude  of  oscilla¬ 
tions  and  maximum  tension  of  the  flexible  coupling  over  the  period 

the  system  described  by  equations  (3)  and  (4)  are  given  in 
Fig.  2  and  3  by  solid  lines.  The  solution  of  equation  (3)  is  described  by  the 
Adams  method  with  zero  initial  conditions. 


Fig.  2.  Amplitude -frequency  characteristics  of  vibra¬ 
tion  protection  systems:  1-R=*=1,  2-12,  3- 
N  -  0.1,  4  -  10. 


As  is  evident  from  Fig.  2,  with  increasing  M  the  amplitude  of  constrained 
vibrations  of  the  stabilized  body  decreases.  Growth  in  the  parameter  M  may  be 
achieved  by  increasing  the  attached  mass  of  liquid  that  is  placed  in  motion  by 
the  body  being  stabilized.  With  increase  in  the  ratio  R  -  R1/R2,  the  ampli¬ 
tude  of  the  resonance  swings  diminishes.  Increase  in  parameter  R  may  be 
achieved  by  improving  the  flow  around  the  auxiliary  body.  Introduction  of  an 
elastic  coupling  ensures  that  the  system  is  placed  in  the  specified  position 
of  static  equilibrium  with  no  perturbing  force.  Lessening  the  stiffness  of 
the  elastic  coupling  reduces  the  lower  limit  of  the  working  frequency  range  of 
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vibration  protection.  The  influence  of  the  nonlinearity  of  the  damping  on  the 
dynamics  of  the  stabilized  body  is  manifested  by  appearance  of  a  second 
resonance  swing,  the  amplitude  of  which  is  much  less  than  that  of  the  first 
resonance  swing  and  diminishes  as  the  parameters  R  and  M  increase.  The 
increase  in  amplitude  of  the  tension  T  in  the  flexible  string  with  rise  in  the 
relative  frequency  q  (Fig.  3)  limits  the  upper  limit  of  the  working  frequency 
region  in  accordance  with  the  condition  T  >  0  (the  flexible  coupling  only 
works  under  tension).  Even  so,  however,  the  working  frequency  range  of  this 
vibration  protection  system  is  broader  than  that  of  the  familiar  ones  [3,  4]. 


The  case  when  the  compensator  (auxiliary  body)  is  located  in  air  (R1/R2  -*  ®) 
is  described  by  equations : 


i+N  i  s  1 5/  +|/g'=  [  l/g’-2/  (ilf+1)  ]sin  T, 

o  -  —  =  1  +  [m,'hQ^i+R,h^Q%  I  i  i  t)  ], 

rriig  m,  g 


Tlie  results  of  a  numerical  analysis  of  equation  (5)  are  presented  in  Fig.  2  by 
the  dot  and  dash  lines.  The  curves  of  the  amplitude  of  vibrations  and  the 
relative  tension  a*  against  the  frequency  q  are  similar. to  the  respective 
curves  for  the  preceding  case .  We  notice  that  the  dynamic  process  becomes 
aperiodic  at  large  values  of  parameters  M  and  N. 

Let  us  consider  a  low-frequency  decoupling  layout  (Fig.  4)  in  which  the 
inertial  force  of  rotation  masses  is  used  as  compensation.  The  stabilized 
body  of  mass  m^^  is  placed  in  liquid  and  joined  by  means  of  a  flexible 
nonstretching  string  to  the  piston  of  a  hydraulic  cylinder  with  cross 
sectional  area  Sq,  arranged  on  a  movable  platform,  such  as  a  vessel.  The 
cavities  of  the  cylinder  above  and  below  the  piston  communicate  respectively 
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with  the  central  and  peripheral  cavities  of  a  hydraulic  cylinder  of  cross 
section  area  S2,  rotating  at  velocity  co,  in  which  are  placed  the 
"compensating"  masses  m2.  Both  cylinders  and  the  lines  which  join  them  are 
fillod  with  an  incompressible  fluid.  The  velocity  of  rotation  is  chosen  to 
ensure  static  equilibrium  of  the  system  when  the  compensators  are  placed  in 
the  middle  of  their  working  stroke  in  the  hydraulic  cylinder.  The  piston  of 
area  Sq  is  connected  to  the  hydraulic  cylinder  by  means  of  elastic  elements. 
Self  regulation  of  the  length  of  the  coupling  between  the  object  and  the 
source  of  perturbations  is  achieved  by  changing  the  position  of  the  piston  of 
area  Sq  relative  to  the  hydraulic  cylinder  in  which  it  is  mounted.  Rotation 
of  the  hydraulic  cylinder  with  the  "compensating"  masses  is  accomplished  by  an 
external  energy  source. 


Fig.  4.1-  stabilized  body,  2  -  piston, 
3  -  nonmoving  cylinder,  4  -  rotating 
cylinder,  5  -  "compensating"  masses,  6  - 
hydraulic  lines,  7  -  elastic  elements. 


The  equation  of  force  balance  for  the  stabilized  body  is  written  as: 

*^0+^ {XQ-Xi) ,  (5) 

where  pj^  and  P2  are  the  pressure  of  fluid  in  the  hydraulic  cylinder  below  and 
above  the  piston  of  area  Sq,  respectively. 

The  equation  of  force  balance  for  the  "compensating"  bodies  is: 

iP>—Pi)S2=mi’Af-(i)*mi’  (r+Ar)  +f^  ( Ar) ,  ^  ^  ^ 

where  m2  is  the  effective  mass  of  the  compensating  bodies,  taking  into 
account  the  attached  mass  of  fluid;  P3,  p^  is  the  pressure  of  the  fluid  in  the 
rotating  cylinder  in  its  central  and  peripheral  cavities ,  respectively;  r ,  Ar 
is  the  static  radius  of  rotation  of  the  bodies  and  the  dynamic  increment  in 
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the  radius;  fi(Ar)  is  a  dissipative  function,  characterizing  the  resistance  to 
displacement  of  the  bodies  in  the  hydraulic  cylinder. 

The  flow  of  fluid  in  the  hydraulic  lines  joining  the  peripheral  cavities  of 
the  rotating  hydraulic  cylinder  to  the  cavity  of  the  nonmoving  hydraulic 
cylinder  beneath  the  piston  and  the  central  cavity  of  the  rotating  hydraulic 
cylinder  to  the  cavity  of  the  nonmoving  hydraulic  cylinder  above  the  piston  is 
described  by  the  Bernoulli  equations  (not  considering  the  height  at  which  the 
cross  sections  are  located) : 


P3  1  (Ar)'_  Pi 

,  {Xa-Xi)-  ,  j. 

2g 

Pi  j  _ 

p,  {^fY^ 

V  2^  ’ 

(8) 


where  7  is  the  specific  gravity  of  the  fluid;  F^,  F2  are  dissipative  functions 
describing  the  decrease  in  pressure  head  in  the  lines. 

The  equation  of  continuity  for  the  fluid  in  the  lines  is: 


^rS2=^(Xi—Xi)So. 

From  (6)  using  (7)>(9),  we  get  the  equation: 


(9) 


\xi  |xl+C;,l=A:v7^2*  (f.  —.f,)  — 

-~k~  (Xo—Xi) 


(10) 


where  k  =  S0/S2. 

The  inertial  and  dissipative  forces  of  the  fluid  may  be  disregarded,  since 
they  are  small  in  comparison  to  the  force  of  hydrodynamic  resistance  of  the 
body  being  stabilized,  i.e., 

/roW  ix,-z,)  ^  ^  F ^  ^ 


Using  definitions  (2)  and  noting  (11),  we  shall  write  the  equation  of 
vibrations  (10)  and  the  expression  for  the  tension  of  the  string  in  the 
nondimens ional  form: 


i+N  1111/  (M+k^)  +|/g^=  [  i/q^-kV  (M+F)  ]  sin  x, 


T 

m,  g 


g|i+C/i(g-siii  T)]/m,'g+l. 


(12) 
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The  results  of  numerical  analysis  of  equation  (12)  are  given  in  Fig.  2  by 
dotted  lines . 

Thus,  the  dynamics  of  the  above  vibration  protection  systems  is  described  by  a 
single  class  of  equations,  and  the  results  of  their  analysis  are  comparable. 
The  investigations  demonstrate  the  possibility  of  ensuring  effective  low- 
frequency  decoupling  of  an  object  from  the  low-frequency  perturbations  caused, 
for  example,  by  the  rolling  of  the  platform  of  the  ocean  exploration 
facilities . 
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the  first  three  paragraphs.] 

[Text] 

One  of  the  drawbacks  of  internal  combustion  automobile  engines  is  that  the  gas 
distribution  phases  do  not  correspond  to  their  optimal  values.  To  eliminate  this,  it 
is  desirable  to  regulate  the  phases  of  gas  distribution  in  the  process  of  engine 
operation. 

Most  promising  is  an  additional  device  before  the  intake  valve. 

Based  on  calculations,  the  requirements  have  been  formulated  for  a  regulating 
cutout  element  before  the  intake  valve. 


One  of  the  drawbacks  of  exiting  automobile  engines  is  the  lack  of  correspondence 
between  the  phases  of  gas  distribution,  especially  the  end  intake  phase,  and  optima! 
values.  This  holds  true  in  the  majority  of  speed  and  load  modes.  By  shortening 
the  phase  one  can  significantly  expand  those  areas  of  speed  characteristics  with  a 
high  torque  and  accommodation  coefficient.  But  at  high  frequencies  of  crankshaft 
rotation  the  torque  and  power  decrease  in  comparison  to  the  torque  and  power 
reached  in  a  longer  end  intake  phase. 

So  began  the  attempt  to  change  the  phases  of  gas  distribution  during  engine 
operation.  However,  only  the  end  intake  phase  and  the  begin  exhaust  phase  could 
be  varied  over  a  wide  range,  since  their  times  are  not  limited  by  the  meeting  of  the 
piston  with  the  open  valve. 

In  recent  years,  as  progress  in  engine  characteristics  has  accelerated,  simple  designs 
have  been  exhausted,  and  new  designs  found  which  were  suitable  for  mass 
production,  it  became  possible  for  industry  to  use  devices  to  implement  near  optimal 
phases  of  gas  distribution  during  engine  operation  in  each  speed  mode  (.Alpha-Romeo, 
Nissan,  Honda,  Mercedes). 

Devices  to  change  the  phases  of  gas  distribution  can  be  divided  into  several  groups 
according  to  their  characteristics.  These  are  devices  for  stepped  and  stepless 
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change  in  the  angular  position  of  the  camshaft  relative  to  the  crankshaft  with 
mechanical,  hydraulic,  electrical,  or  combined  control  [1], 

The  simplest  of  the  known  devices  is  a  stepless  means  of  changing  the  angular 
position  of  the  camshaft  relative  to  the  crankshaft.  Screw,  chain,  geared  drive,  or 
differential  mechanisms  or  devices  are  used  for  this.  ’ 

The  simple  design  of  these  devices  conflicts  with  their  principal  drawback, 
simultaneous  change  in  the  phase  of  the  beginning  and  end  of  intake.  Practical 
means  of^varying  the  phases  in  this  way  do  not  exceed  30-40”  of  the  rotation  of  the 
>.rankshaft  because  one  must  prevent  the  collision  of  open  valves  and  pistons.  But 
even  such  a  relatively  small  change  in  phases  yields  positive  results:  one  can 
increase  the  filling  of  cylinders  at  low  rotation  frequencies  and  increase  the  supply 
of  torque  and  the  power  per  liter. 

A  more  complex  means  of  stepless  regulation  of  gas  distribution  phases  is  based  on 
the  use  of  a  camshaft  consisting  of  two  mutually  joined  parts.  In  each  part  (the 
internal  and  external)  are  cams  whose  sides  touch.  The  mechanism  controlling  the 

mutual  movement  of  the  cams  makes  it  possible  to  separate  them  by  an  angle  of  up 
to  44  .  ay 


Another  design  was  developed  by  Fiat.  Fiat  created  an  mechanism  with  cams  which 
were  wider  than  the  surface  of  the  push  rod  sliding  along  it  and  with  a  variable 
profi.e  along  the  axis.  The  regulation  was  done  by  shifting  the  camshaft  along  the 
axis.  In  a  means  of  stepped  phase  regulation  patented  by  the  Matsuda  firm,  the 
camshaft  has  two  sets  of  cams  with  different  profiles  (for  low  and  high  rotation 
irequencies/.  Each  valve  of  the  mechanism  is  comprised  of  two  rockers,  one  of 
which  has  a  constant  kinematic  link  with  the  valve.  The  other  comes  into  contact 
by  blocking  in  the  transition  to  a  high  speed  mode. 

In  a  separate  group  are  additional  cutout  devices  in  the  intake  tract  which,  in 
addition  to  regulating  phase,  decrea.se  the  throttle  losses  which  are  inherent  in 
engines  with  quantitative  regulation  and  a  rotation  which  increases  sharply  as 
frequency  decreases.  Volkswagen  has  used  a  rotating  slide  valve  (Figure  1)  before 
tile  intake  valve  in  an  experimental  engine. 


Figure  1.  1.  intake  channel;  2.  rotating  slide  valve;  3.  intake  valve. 

The  slide  valve,  which  is  driven  by  a  device  which  varies  its  angular  position 
relative  to  the  crankshaft,  may  cut  off  the  intake  channel  at  any  moment  (without 
limit)  ana  cnange  the  effective  phase  of  the  end  of  intake  by  shortening  it  in 
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relation  to  the  corresponding  phase  of  the  main  intake  valve.  Rut  since  the  slide 
valve  rotates  with  a  constant  angular  speed  equal  to  the  speed  of  the  main 
camshaft,  its  window  does  not  have  fixed  open  and  closed  positions,  which  inevitably 
leads  to  an  increase  in  throttle  losses, 

VToreover,  the  slide  valve  has  a  long  developed  cylindrical  surface.  If  one  is  to 
provide  a  small  clearance  in  the  slide  valve  in  relation  to  the  immobile  parts  in 
contact  with  it,  this  leads  to  the  passage  of  a.  live  charge  through  the  contact  while 
"advancing"  the  cutout  of  the  cylinder  an  the  intake  manifold. 

The  use  of  such  a  device  reveals  great  possibilities  for  controlling  the  actual  phase 
of  the  end  of  intake,  from  delaying  it  to  significantly  advancing  it  in  relation  to 
lower  dead  center.  In  the  latter  case,  one  can  sharply  decrease  the  charge  of  the 
cylinders  without  throttling  in  the  carburetor,  that  is,  without  creating  a  significant 
vacuum  in  the  intake  collector.  Indeed,  the  prerequisites  for  theoretically  throttle¬ 
less  quantitative  regulation  are  created.  However,  practically,  due  to  aerodynamic 
imperfection  of  the  “broken'’  line  of  the  intake  channel  from  the  carburetor  to  the 
vaive  due  to  partial  mismatching  of  the  axes  of  circumference  in  the  slide  valve 
and  in  the  nozzle,  the  throttle  losses  are  partially  preserved. 

An  interesting  feature  of  the  ideal  process  of  filling  the  cylinder  in  the  presence  of 
this  type  of  device  is  the  near  adiabatic  decrease  in  the  charge  temperature  when 
its  entrance  to  the  cylinder  is  closed  early  (in  relation  to  lower  dead  center).  The 
closer  the  additional  cutout  device  to  the  valve,  the  higher  the  degree  of  cooling. 
A  calculation  using  a  very  simple  expression  of  the  polytrope  when  its  index  is 
close  to  1,4  shows  that  when  the  entrance  of  charge  into  the  cylinder  is  closed  40*" 
before  lower  dead  center,  its  temperature  decreases  by  30",  and  when  it  is  closed  45" 
before  lower  dead  center,  temperature  decreases  by  50".  Such  substantial  cooling 
aids  an  intense  decrease  in  the  tendency  to  detonate,  its  prevention  despite  the 
very  early  onset  of  ignition  and  the  low  octane  number  of  the  fuel.  Thus,  one  can 
sharply  reduce  the  requirement  on  the  antidetonation  properties  of  benzine  without 
worsening  fuel  economy  in  the  most  characteristic  traveling  modes.  However,  the 
design  proposed  by  Volkswagen  is  far  from  ideal. 

I.et  us  examine  the  effect  of  deviations  from  the  theoretical  cycle  on  energy  losses. 

1.  The  presence  of  a  gap  in  the  contact  with  the  cylinder  head  leads  to  the 
entrance  of  fuel  into  the  cylinder  when  the  slide  valve  is  closed,  but  when  the 
intake  valve  is  still  open.  Considering  the  value  of  the  entrance  of  fuel  in  the 
cylinder  through  the  gap,  we  get  its  dependence  on  the  value  of  the  gap  (Figure  2, 
curves  1  and  2).  It  is  clear  from  the  graph  that  even  a  small  increase  in  the  gap 
makes  the  course  of  the  working  cycle  much  worse.  As  the  frequency  of  rotation 
of  the  engine  increases,  this  correlation  increases.  But  due  to  the  great  length  of 
the  slide  valve  it  is  rather  difficult  to  provide  a  small  gap  in  its  contact  with  the 
casing. 

2,  The  presence  of  a  “parasite”  volume  between  the  slide  valve  and  the  intake 
valve  leads  to  a  decrease  in  the  degree  of  regulation  of  the  engine,  which  is  linked 
with  the  difference  in  the  volumes  of  polytropic  expansion  in  the  cycle  of  intake 
and  subsequent  compression.  By  using  the  equation  of  the  polytrope  we  obtain  the 
dependence  of  the  theoretical  losses  of  this  type  on  the  “parasite”  volume  (Figure 
3).  It  follows  from  these  relations  that  the  “parasite”  volume  must  be  decreased  to 
about  10%  of  the  cylinder  volume. 
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Figure  1.  radial  gap  of  O.OS  mm;  2,  radial  gap  of  0.1  mm;  closure 
equal  to  the  time  for  the  camshaft  to  rotate  (3)  10”,  (4)  a 

Q'/Q  min 


time  of  the  gate 
Ogap/G^yt,  %;  b. 


Figure  3.  “Parasite”  volume  composed  of  the  volume  of  the  cylinder;  1,  10;  2,  20;  3, 
30;  4  engine  without  an  additional  device  at  the  intake.  a.  P,  kg/cm^  *  b! 

^rot.  ira?iksh-.x/t' 

3.  The  final  rate  of  cutout  of  the- intake  line  by  the  slide  valve  is  proportional  to 
the  .speed  of  the  crankshaft.  The  energy  losses  are  due  to  throttling  of  the  flow  in 
cross  sections  defined  by  the  coverage  of  the  valves  in  the  slide  valve,  the  “break” 
in  the  axial  line  of  the  intake  tract,  and  the  presence  of  breaks  in  the  flow  when  it 
flows  around  sharp  edges. 


36 


To  estimate  the  loss,  Pinite-difrerence  expressions  were  used  which  were  based  on 
experimental  and  theoretical  studies.  In  relation  to  the  Reynolds  number  [2] 


APr  = 


A_ 

Re 


1  -  g|  e.,. 


p 

^  i 


where 


A  -  120 


''l  4-  0.5D^  (1  -t-  sin  d~)1 
(1  -  Dl  sin  ! 


Re  is  the  Reynolds  number,  is  the  given  diameter  of  the  gate,  S  is  the  angle  of 
rotation  of  the  gate,  is  a  coefficient  which  considers  the  limit  losses,  p  is  the 
density  of  the  medium,  Q  is  the  volume  expenditure,  and  F  is  the  cross  sectional 
area. 

Calculations  have  made  it  possible  to  establish  the  dependence  of  the  throttle  losses 
on  the  rate  of  closure  of  the  cutout  element  (Figure  2,  curves  3  and  4).  Thus,  the 
effect  of  the  rate  of  the  cutout  of  the  intake  line  on  filling  is  significant,  from  4% 
to  10%. 

f.ong  before  Volkswagen,  the  Scientific  Research  Automobile  and  Auto  Motor 
Institute  received  an  inventor’s  certificate  on  a  means  of  using  an  effective  degree 
of  compression  without  the  drawbacks  which  are  characteristic  of  devices  with  a 
slide  valve  cutout  mechanism  [3].  In  comparison  to  the  latter,  the  device  proposed 
by  the  Institute  sharply  reduces  the  work  of  the  pumping  strokes,  decreases 
throttle  losses,  reduces  flow  through  the  closed  cutoff  device,  reduces  the  volume 
between  it  and  the  main  valve,  and  reduces  the  energy  expenditures  for  the  drive. 

Conclusions 

The  use  of  an  additional  cutout  device  in  the  intake  tract  before  the  main  intake 
valve  for  quantitative  regulation  of  an  internal  combustion  engine  requires  the 
following  conditions:  1)  a  small  gap  in  its  contact  with  the  casing;  2)  the  volume 
between  the  additional  device  and  the  intake  valve  should  not  exceed  10%  of  the 
cylinder  volume;  3)  the  rate  of  closure  of  the  cutout  device  should  be  shorter  than 
the  time  for  the  crankshaft  to  rotate  40\ 
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ABSTRACT 

This  article  describes  the  designs  and  principles  of  action  of  sensors  and  primary 
transformers  for  the  measurement  of  the  strength  of  friction  on  a  wall  in  reverse 
flows  with  the  direction  of  movement  varying  over  time.  The  work  of  the  examined 
measurement  devices  is  based  on  the  use  of  the  thermoanemometer  method  to 
determine  the  instantaneous  strength  of  friction  and  the  use  of  the  heat  trace  from 
the  sensitive  element  of  the  thermoanemometer  to  determine  the  instantaneous 
strength  of  the  action  of  the  force  of  friction.  Some  results  are  presented  of  the 
testing  of  the  method  in  a  stalled  flow  before  a  ledge  in  a  channel. 


Complex  turbulent  flows  with  stalled  zones  are  encountered  in  the  flow  tract  of 
energy  installations.  A  characteristic  of  these  flows  is  the  nonstationary  position 
of  the  stall  and  junction  points.  Another  characteristic  is  random,  instantaneous 
reversal  in  the  direction  of  movement  of  the  flow  in  the  vicinity  of  these  points. 
Common  methods  of  measuring  the  strength  of  friction  in  the  wall  using  a  floating 
element,  Preston  tubes,  Stanton  “deflectors”,  and  Gardon  type  heat  sensors  [1]  are 
not  suited  to  the  investigation  of  such  flows.  A  pulse  anemometer  sensor  is  known 
[2]  which  can  measure  the  average  value  of  the  strength  of  friction  in  a  reverse 
flow.  However,  it  has  a  very  low  frequency  range  and  cannot  measure 
instantaneous  pulsations  in  friction,  the  frequency  of  the  change  in  the  direction  of 
movement,  and  other  characteristics. 

Below,  two  thermoanemometric  methods  are  described  for  the  measurement  of  the 
strength  of  friction  in  reverse  flows.  In  the  first  method,  a  sensor  is  used  with 
three  parallel  thin  threads  located  in  one  plane  at  a  distance  H  =  0.05  0.3  mm 
from  the  surface  of  the  wall.  The  central  thread  1  (Figure  1)  of  the  sensor  is  the 
sensitive  element  of  the  thermoanemometer,  and  determines  the  instantaneous  value 
of  the  strength  of  friction  at  the  wall 

T  =  fji{du/dy)y.Q  «  (J.Uh/H,  (1) 
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Figure  1.  Schematic  of  the  three-thread  sensor  and  primary  transformer,  a.  out. 

where  //  is  the  dynamic  viscosity,  and  u  is  the  speed  of  the  flow  at  distance  ff  from 
the  wall.  Two  other  lateral  threads,  2  and  .3,  of  the  sensor  are  thermometers  of 
resistance.  Their  purpose  is  to  determine  the  instantaneous  direction  of  movement 
of  the  flow,  or  the  instantaneous  direction  of  action  of  the  force  of  friction  along 
the  heat  trace  from  the  heated  thread  of  the  thermoanemometer.  These  threads  are 
2.5-5.0  /y.m  in  diameter  in  fine  holders  0.2  mm  in  diameter  4  located  at  0.2-0.4  mm 
from  the  central  thread  on  both  sides. 

Tn  the  first  method,  a  standard  thermoanemometer  5,  linearizer  6  and  correcting 
adapter  are  used  as  primary  transformers  to  measure  friction.  The  adapter  consists 
of  a  differential  amplifier  7,  comparator  8,  electronic  switch  9,  and  an  inverter- 
repeater  10.  The  input  of  the  differential  amplifier  is  connected  with  the  diagonal 
of  the  resistive  bridge.  Resistance  thermometers  2  and  3  are  connected  at  the 
bridge’s  shoulder. 

When  thread  1  is  not  heated,  the  bridge  is  balanced.  When  the  thread  is  heated,  its 
heat  trace,  depending  on  the  direction  of  the  flow,  falls  on  one  of  the  thermal 
resistances.  This  creates  either  a  positive  or  negative  voltage  which  unbalances 
the  bridge.  This  voltage  is  amplified  by  the  differential  amplifier  and  enters  the 
comparator,  which  has  a  zero  response  threshold.  The  comparator  forms  signals  to 
control  the  electronic  switch  of  the  inverter-repeater.  At  the  input  of  the 
inverter-repeater  the  linearized  output  signal  of  the  thermoanemometer  enters. 
This  signal  is  proportional  to  the  instantaneous  value  of  the  strength  of  friction  at 
the  wall.  Depending  on  the  direction  of  the  flow  the  inverter-repeater  either 
inverts  or  repeats  the  phase  of  the  input  signal.  Thus,  a  correlation  is  established 
between  the  output  signal  of  the  thermoanemometer  and  the  instantaneous  value  of 
the  strength  of  friction  at  the  wall  considering  the  direction  of  movement  of  the 
flow. 

The  other  method  uses  a  two-thread  sensor.  Rach  of  the  threads  is  the  sensitive 
element  of  the  thermoanemometer.  The  structure  and  size  of  the  two-thread  sensor 
is  approximately  the  same  as  the  three-thread  model.  Threads  1  and  2  (Figure  2) 
are  parallel  and  0.1-0.4  mm  apart  from  each  other.  The  primary  transformer  is  two 
thermoanemometers  3  and  4  with  linearization  units  5  and  6  and  a  correcting  adapter. 
The  adapter  consists  of  a  differential  amplifier  7,  a  matching  cascade  8,  an  inverter 
9,  two  electronic  switches  10  and  11,  a  comparator  12,  and  summer  13.  The 
instantaneous  direction  of  movement  of  the  flow  is  determined  from  the  difference 
of  the  readings  of  thermoanemometers  3  and  4.  The  thermoanemometer  in  which  the 
sensitive  element-thread  is  in  the-  heat  trace  reduces  the  output  signal.  The 
instantaneous  strength  of  friction  is  determined  by  the  same  thermoanemometer 
whose  sensitive  element  is  not  in  the  heat  trace  or  is  first  in  the  flow. 

Let  us  briefly  describe  the  working  principles  of  the  correcting  adapter  (Figure  2). 
The  linearized  signals  of  thermoanemometers  3  and  4  enter  differential  amplifier  7, 
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Figure  2.  Schematic  of  the  two-thread  .sensor  and  primary  transformer. 

which  compares  the  input  signals  and  controls  comparator  12.  The  comparator  forms 
signals  for  the  control  of  switches  10  and  11.  Depending  on  the  direction  of 
movement  of  the  flow,  switch  10  is  opened  and  the  signal  of  thermoanemometer  ."1 
passes  to  the  output  of  the  device,  or  switch  11  is  opened  and  the  inverted  signal 
of  thermoanemometer  4  passes  to  the  output. 

The  electrical  schematics  of  the  correcting  adapters  which  determine  the 
instantaneous  direction  of  movement  of  the  flow  in  the  two  methods  examined  are 
presented  in  [3].  We  implemented  these  adapters  in  the  form  of  instruments  and 
tested  them  together  with  sensors  for  their  working  capacity.  The  standard 
thermoanemometers  used  were  equipment  of  the  DISA-55M  system,  which  made  it 
possible  to  make  measurements  of  instantaneous  flow  speed  with  a  pulsation 
frequency  of  up  to  50  kHz. 

The  testing  showed  that  the  three-thread  sensor  makes  it  possible  to  record  with 
high  accuracy  changes  in  the  direction  of  movement  of  the  flow  in  the  opposite 
direction  with  a  frequency  of  up  to  100  Hz.  The  lateral  threads  have  a  slight 
effect  on  the  reading  of  the  central  thread  of  the  thermoanemometer.  The 
probability  of  a  possible  failure  of  the  thermal  trace  to  fall  on  the  lateral  thread 
for  the  de-scribed  sensor  construction,  even  directly  in  the  vicinities  of  the  stall 

and  junction  points  of  the  flow,  did  not  exceed  1%.  which  has  virtually  no  effect  on 
the  results  of  the  measurements.  The  testing  or  the  two- thread  sensor  showed  a 
larger  range  of  frequencies  of  flow  reversal,  up  to  500  Hz.  This  increase  in  the 
frequency  range,  in  comparison  to  the  three-thread  sensor,  is  explained  by  the  fact 
that  in  that  method  both  threads  operate  in  constant  temperature  mode.  The 
studies  of  many  types  of  subsonic  stalled  flows  that  we  performed  showed  that  the 
frequency  of  reversal  of  flow  does  not  exceed  40-70  Hz. 

Testing  of  the  described  sensor  structures  also  showed  that  the  sensors  must  be 
friction  graded.  There  are  a  number  of  problems  associated  with  determining 
friction  directly  from  formula  (1).  First,  in  manufacturing  of  the  sensor  it  was 
difficult  to  maintain  the  identical  distance  II  along  the  entire  thread  length. 
Second,  one  had  to  consider  the  effect  of  the  closeness  of  the  wall  on  the  reading 
of  the  thermoanemometer.  Another  means  of  determining  the  strength  of  friction  at 
the  wall  is  with  direct  friction  grading  in  a  common  boundary  layer  on  film,  which 
is,  in  our  opinion,  less  labor-intensive  and  more  reliable. 

It  should  also  be  noted  that  in  friction  sensors,  instead  of  threads  one  may  use 
metal  films  applied  to  a  special  substrate.  The  u.se  of  films  makes  it  pos.sib1e  to 
exclude  the  disturbing  effect  of  the  sensitive  elements  of  the  sensor  on  the  flow. 

T.et  us  present  several  results  of  the  study  of  flow  in  the  region  of  a  flow  stall 
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before  a  ledge.  These  results  were  obtained  using  a  three-thread  friction  sensor. 
The  experiments  were  done  in  a  subsonic  aerodynamic  tube  which  operated  in 
exhaust  mode.  The  working  part  of  the  tube  was  a  rectangular  channel  140  mm  high 
and  360  mm  wide.  The  ledge  was  42  mm  high  and  was  located  across  the  full  width 
of  the  lower  wall  of  the  channel.  The  speed  of  flow  before  the  ledge  was  regulated 
from  5  to  36  m/s.  The  collection  and  processing  of  information  during  the 

experiments  was  done  using  a  measurement-computing  complex  based  on  a  SM-4 

computer  using  the  Turbulence  automated  measurement  system  [41. 

As  a  result  of  the  processing  of  experimental  results  data  were  obtained  on  the 
average  and  pulsation  components  of  the  strength  of  friction  r,  the  alternation  7, 
and  the  frequency  P  of  change  in  the  direction  of  movement  of  the  flow  over  time. 
The  alternation  was  determined  as  a  portion  of  time,  during  which  the  direction  of 
movement  of  the  medium  at  the  mea.surement  point  coincided  with  the  direction  of 
movement  of  the  external  flow. 

Figure  3  shows  the  schematic  of  flow  and  some  results  of  the  measurements  of 
characteristics  of  flow  in  the  stall  region  before  the  ledge.  Analysis  of  the  results 
of  measurements  showed  that  two  flow  modes  were  possible  before  the  ledge:  with 
the  generation  of  a  small  recirculation  zone  before  the  stall  region  and  without  it. 
In  the  first  case  the  friction  at  the  wall  changes  sign  four  times,  and  in  the  second, 

two  times.  The  maximum  value  of  the  average  strength  of  friction  at  the  wall  is 
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Figure  3.  Structure  of  the  flow  in  the  region  of  flow  stall  before  the  ledge.  Solid 
line,  =  15  m/s,  ^  0.013;  »  0.116.  Dashed  line,  u-  = 

15  m/s. 
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obsserved  approximately  at  the  midpoint  between  the  stall  point  and  the  base  of  the 
ledge,  and  the  maximum  pulsations  in  friction  in  the  region  of  the  flow  stall 

from  the  wall,  where  f  =  0.  In  the  regions  of  stall  and  junction  of  the  flow,  where 
the  average  strength  of  friction  at  the  wall  was  equal  to  zero,  7  50%,  and  F 

reached  maximum  values.  This  result  indicates  that  the  points  of  stall  and  junction 
of  the  flow  may  be  determined  from  the  alternation  values  and  the  frequency  of 
change  in  direction  of  the  movement  of  the  flow.  Friction  grading  of  the  .sensors  is 
not  necessary  to  measure  alternation  and  frequency. 

Thus,  the  suggested  methods  significantly  expand  the  region  of  use  of 
thermoanemometry  to  measure  the  strength  of  friction  in  reverse  flows. 
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